Very High Mach Number Electrostatic Shocks in Collisionless Plasmas 
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The kinetic theory of collisionless electrostatic shocks resulting from the collision of plasma slabs 
with different temperatures and densities is presented. The theoretical results are confirmed by self- 
consistent particle-in-cell simulations, revealing the formation and stable propagation of electrostatic 
shocks with very high Mach numbers (M 3> 10), well above the predictions of the classical theories 
for electrostatic shocks. 
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The collision of clouds of plasma with different prop- 
erties (temperature, density, composition etc.) is a sce- 
nario quite common in nature. For instance, during su- 
pernovae explosions, large quantities (10 solar masses) of 
high temperature plasma ( T ~ 10 6-8 K) are ejected into 
the interstellar medium (n ~ 1 cm" 3 , T ~ 10 2 " 4 K) [HQ, 
and plasma cloud collisions are at the core of the fireball 
model for gamma ray bursters Q. Plasma cloud colli- 
sions also occur when the solar wind interacts with the 
Earth Magnetosphere, or when it encounters the inter- 
stellar medium in the heliosphere region 0]. In the lab- 
oratory, such scenarios appear during the laser induced 
compression of plasma foils in solid targets []| . 

The collision of plasma shells leads to the onset of 
plasma instabilities and to the development of nonlinear 
structures, such as solitons, shocks and double layers 0. 
In the absence of an ambient magnetic field, the shock 
waves are electrostatic HI , and the dissipation is pro- 
vided by the population of electrons trapped beyond the 
shock front [fjlij and, for stronger shocks, by the ion re- 
flection from the shock front [lOj. Whilst the properties 
of shocks induced by collision of identical plasma shells, 
or by compression of plasma clouds, have been exten- 
sively studied in the past 0, IE 0, S H E3, Ell E3, E3, d , 
the properties of the electrostatic shock waves formed 
during the collision of diverse plasma slabs of arbitrary 
temperature and density are rather unexplored |l5| . The 
theory for electrostatic shocks induced by impact of iden- 
tical plasma shells predicts an absolute maximum Mach 
number M max ~ 3 (or, when ion reflection and thermal 
effects are included, M^ ax ~ 6). However, collisionless 
shock waves with Mach numbers ranging between 10 and 
10 3 have been observed in many astrophysical scenarios 
, and very large Mach number cosmic shock waves are 
thought to play a crucial role in the evolution of the large 
scale structure of the Universe 0, 0] . 

In this Letter, we present a kinetic theory describing 
the properties of the very high Mach number (M ^> 10) 
laminar shock waves arising from the collision of slabs 



of plasma with different properties (temperatures, den- 
sities) , and in the absence of an ambient magnetic field. 
We demonstrate that the shock properties are strongly 
influenced by Q, the ratio of the electron temperatures 
in the two slabs, and by T, the ratio of the electron den- 
sities in the two slabs. The analysis shows that when the 
electron temperature T c r of the downstream slab (R) is 
higher than the electron temperature T e l of the upstream 
slab (L), the shock waves can have very large Mach num- 
bers, which are otherwise not supported by isothermal 
plasmas 0, 0]. The model predicts that the maxi- 
mum allowed Mach number increases with 0, without 
an absolute upper limit. The theoretical results are con- 
firmed by one-dimensional (ID) self-consistent particle- 
in-cell simulations, demonstrating the formation and the 
stable propagation of electrostatic shock waves with very 
large Mach numbers (M ~ 20). 

The shock transition region is modeled as a planar one 
dimensional double layer, which is stationary in the ref- 
erence frame of the shock; the electrostatic potential in- 
creases monotonically from (j) — at x — xl to <p = cf>o 
at x — xr, as shown in Figure ^ The one dimensional 
treatment should hold for propagation distances smaller 
than the transverse dimension of the shock. The struc- 
ture of the double layer is maintained by two popula- 
tions of free electrons (from the L and R slabs), a pop- 
ulation of trapped electrons (from the R slab), and a 
population of cold ions. The model considers the con- 
tribution of the trapped particles self-consistently, by 
treating the electrons in a kinetic fashion 0, 0- To 
maintain a steady state, the number of electrons in unit 
time dt, with velocity between [v e ,v e + dv e ] and po- 
sition between [x,x + dx], must balance the electrons 
injected from the left slab (L), at i = xl, with ve- 
locity between [vt,,vl + dvr,], and the electrons from 
the right slab (R), at x = xr, with velocity between 
[vr> vr + dvn]. We can express the electron distribution 
function f e (x,v e ), at any given point x, as a function of 
the electrons injected at the left and right boundaries, 
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which follow known distribution functions. The electron 
velocity, which results from conservation of energy, can 
be written as v e = + ^ = -y^ + MizM de- 
pending if the particles are arriving from the left bound- 
ary or from the right boundary; here m e is the electron 
mass, and e is the elementary charge. Assuming that 
only electrons with positive velocity enter from the slab 
L, and only electrons with negative velocity enter from 
the slab R, we obtain: 



f e (x,v e )dv e 



fh(v L )dv L 



Vr 



fR(vn) d,V! 



(1) 

where and /r(ur) are the electron distribution 

function in the L and R slabs. We assume that the 
electrons in the L slab follow the drifting Maxwell- 
Boltzmann (MB) distribution function |6( /l(i>l) = 
,-(v L -Vi) /V TeL anc j we consider that the elec- 
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tron thermal velocity Vrea = \/ k R T c a / m e is greater 
than the shock speed V*. Here fee is the Boltzmann con- 
stant, while T ea and No a are the electron temperature 
and density in the slab a. In the R slab, the electrons are 
affected by the potential </>o, and their distribution func- 
tion is composed of a free and a trapped part, /r(vr) = 
fnf{v R ) + fnt(vR)- If their kinetic energy is larger than 



the electrostatic energy (i.e. \v R \ > 



are free, and continuously decelerate while moving to- 
wards the left boundary. Following 0, 0, we assume 
that the free electrons can be modeled by the distribution 



function f R f(v R ) 
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e 4>o 



oR . The elec- 



trons with kinetic energy lower than the electrostatic en- 
ergy (i.e. \v R \ < \v c \) are trapped, and are assumed to fol- 
low the flat top distribution function fm(v R ) — T , NaR ^ - 

VTeRW 2.TT 

This idea of describing an electron gas, composed of free 
and trapped particles, by using a MB distribution func- 
tion with a flat-top, has been widely used in the pas t 
to model stationary ion acoustic structures [1, 0, 0, ^3 > 
and it is called 'maximum-density-trapping' approxima- 
tionHQ. 

The density of electrons along the shock can be cal- 
culated by integrating the electron distribution function 
f e (x,v e ) in velocity space. Following Equation QJ, we 
can separate the electron density along the shock in the 
contributions from the L and R slabs. The density of 
the electrons coming from the slab L can be written as 
nh(tp) = -^j L e v Erfc v /(p, where tp = fe e ^ . In our model, 
the electrons coming from the region L are continuously 
accelerated while moving towards the right boundary, 
and are not reflected or trapped by the electrostatic po- 
tential. The density of electrons from the slab R is ob- 
tained using Equation lfT|l. integrating the term contain- 
ing the free electron distribution, fB.f(v R ), in the range 
[— oo, — v c ], and the term containing the trapped elec- 
tron distribution, /Rt(vfi), in the range [— v c , 0], leading 



to n K (ip) = ^E!eV/ e ErfcVW0+ ^^V^l®' where 
T = N or /N l is the density ratio. 

In the present model, the ions are cold, flowing to- 
wards the shock with velocity Vi , and being continuously 
decelerated by the electrostatic potential. The ion den- 
sity is determined by considering the energy and the 
mass conservation equations |l3j . and can be written 
as m = N / (y/l - 2tp/M 2 y where M = V t /V s is the 

ion acoustic Mach number, V s = y / KBT^Jrrii is the ion 
sound speed, rrii is the ion mass, and Nq is the unper- 
turbed ion density in the slab L. The reflection of ions 
is not included in the present model. Such assumption 
is consistent with a double layer solution maintained by 
a population of trapped electrons 0, . The aim of the 
present Letter is to show that the collision of plasma 
slabs, with appropriate temperature and density ratios, 
leads to shock waves with very large Mach numbers. 
Since the occurrence of ion reflection increases even fur- 
ther the velocity of the shocks, the conclusions of the 
present Letter are to be considered as conservatives. 

In the derivation of the ion and electron densities, we 
used the quantities Nor, ./Vol and No, which can now 
be evaluated by applying boundary conditions proper of 
double layers jy, |2( . Using charge neutrality at x = xl 
and at x = x R , we obtain the conditions No = Nql + Nq^ 
and tir(i/?o) + ^l(vo) = n i{ l Po), thus leading to JVol and 
,/Vor as function of the unperturbed ion density No, of 
the normalized potential (po, and of the Mach number 
M . Since it is clear that the electron distribution func- 
tions are always positive, we must apply the physical 
inequalities N l,N r > 0. 

By combining the ion and electron densities with Pois- 
son's equation, and since the dynamics of the electro- 
static potential is analogous to the motion of a particle 
in a potential well '5, we find that the evolution of the 

electrostatic potential is governed by | ( ) +^(<p) = 
|l8| , where the spatial coordinate, \i ls normalized to 
the electron Debye length Xd = \/kRT jj4:Tre 2 Nn, and 
the nonlinear Sagdeev potential 10] is: 



%) = -{p e ( v ,e ! T)-%,M)}. 



(2) 



Here P e (p,0,T) = PfoG.T) - P(<p = 0,0, T) 
is the electron pressure, and Pi(ip,M) = 

M 2 ^1 — y/l — 2ip/M 2 ^j is the ion pressure, normalized 

to N a k B T cL . The pressure term P(<p, 0, T) = P L + Pr 
includes the contribution of the electrons from the slab 

N 



L, Pl(^) 
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and the contri- 



bution of the free and trapped electrons from the slab 
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R, viz Pr(^) = [e«Erfc/§ + 

The trapping potential in Equation (J2J) was obtained 
by assuming \&(0) = 0. From the analogy with particle 
motion, when the Sagdeev potential J2J) is negative the 
electrostatic potential is driven out of equilibrium and 
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the system supports soliton-like structures. On the other 
hand, the conditions of charge neutrality at the bound- 
aries, equivalent to considering d^(ip)/dip\o yVQ = 0, 
assure that (p grows monotonically from to ip with- 
out oscillating back and forth, while the condition 
^(tfo) = ^(0) = assures that po remains bounded 
without growing indefinitely. 

Examining Equation we can thus conclude that 
the system supports a monotonic double layer solution, 
for a given Mach number, only if the electron pressure 
exceeds the ion pressure along the shock, and if both 
coincide in value and slope at the boundaries x = x^ (tp = 
0) and x — x-^{p = ipo). On the other hand, if the 
electrostatic potential exceeds the critical value ip cr — 
M 2 /2, the ion pressure becomes imaginary, and the wave 
'breaks' @. In order to have a steady state solution, we 
must then impose ^>{tp cr ) > The inequality can 

be written in terms of ion and electron pressures in the 
form P e (M 2 /2,0,T) < P\ = M 2 , which recovers the 
same results obtained by previous authors in the limit of 
— * 1 and T — ► 1 |6J, |]j , and that imposes an upper limit 
to the Mach number of the shock waves created during 
the collision of two plasma slabs with temperature ratio O 
and density ratio T. When M 2 > 1 and M 2 > V©, the 
electron pressure P e (ip,Q,Y) can be properly expanded 
around ip — ip cr , and the expression for the maximum 
Mach number can be written as: 




(3) 



Equation © shows that the collision of two plasma slabs 
can give rise to electrostatic collisionless shocks with very 
large Mach number, provided that the two slabs have the 
appropriate temperature (O) and density (T) ratios. As 
far as we know, no electrostatic shock waves with M > 6 
have been predicted before. This difference arises from 
the the fact that the present model includes the variation 
of electron pressure not only as a function of the electro- 
static potential, but also as function of the temperature 
jump, O, between the electrons in the downstream and 
upstream plasma slabs. When the downstream electrons 
have a temperature larger than the upstream electrons, 
the electron pressure in the shock is reduced (few elec- 
trons are trapped), the maximum electrostatic potential 
ip cr can stabilize to a larger value defined by O, and the 
maximum Mach number M max = \J2p cr increases ac- 
cordingly. Figure shows the theoretical prediction of 
the maximum Mach number as a function of O , for differ- 
ent plasma density ratios T. The solid line represents the 
collision of two plasma slabs with equal density (T = 1), 
recovering the classic limit M ~ 3.1 0, El when 8 = 1. 
The condition for the minimum Mach number is found 
by imposing that the Sagdeev potential is negative at its 
minimum. When = 1 and ./Vol = Aor, the minimum 
Mach number is 1, as in the hydrodynamic limit. 

In order to check the consistency of the theoretical 



predictions, we have performed particle- in-cell simula- 
tions using the fully relativistic massively parallel code 
OSIRIS 2.0 [Hj. The ID simulations are performed in 
the reference frame of the slab L, distance is normal- 
ized to c/uipeL, charge to the electron charge e, mass 
to the electron mass m e , and time to l/u) pe L, where 

1/2 

Wp e L = (47re 2 A L/m e ) . The box length is 120 c/ui pe L, 
with 32768 cells, 50 particles per cell per species (4 
species), and the time step is oj pe L<lt = 3.63 x 10 -3 . The 
simulations start at u> pe L t = 0, with the slab L occupying 
the region xlu ps l/c = [0,80], and the slab R occupying 
the region xLu pe x/c = [80,120]. In the simulations, the 
shock is driven by the slab R, both moving to the left; the 
shock and driver velocities are calculated in the frame of 
the slab L. The simulations cover a wide range of param- 
eters, with the driver velocity varying between Mach 2 
and Mach 40. The electrons in the R and L slabs have 
temperatures T c r = 1 keV and T c l = 10 eV, respectively, 
the temperature ratio is O = 100, and the density ratio 
is T = 3. In such conditions of temperatures and densi- 
ties, depending on the velocity of the driver, the theory 
predicts the formation of shock waves with Mach number 
ranging between 10 and 20 (cf. Figure0). Figure [3] shows 
the comparison between the theoretical and the numer- 
ical electrostatic potential ip, as a function of the Mach 
numbers of the shocks observed in the simulations; such 
shocks showed velocities ranging between 10 and 20 times 
the ion sound speed, in excellent agreement with the the- 
ory. It should be noted that, while no shocks are observed 
for drivers (slab R) moving with velocity above Mach 20, 
very high Mach number shocks (M ~ 10 — 20) are cre- 
ated by drivers moving with velocity between Mach 2 
and Mach 20. In the simulations, the value of the elec- 
trostatic potential has been calculated for well developed 
shock structures, i.e. normalized times much larger than 
u pe Lt = 1000. Figure 0] shows the typical phase space 
plxl of the ions from the L slab, at four different time 
steps. The slab R moves towards the left with M = 15, 
and drives an electrostatic shock moving at M = 16. 
The electrostatic potential predicted by the theory for 
the conditions in our simulation, ipx — 123, is in very 
good agreement with the electrostatic potential obtained 
in the simulations, ipo — 125 ± 5. As the shock struc- 
ture propagates, the ions are picked up and accelerated: 
a small fraction is reflected by the electrostatic potential, 
while most of the ions end up with the same speed of the 
shock wave. We have also performed simulations of sce- 
narios with two identical slabs: in these scenarios, when 
colliding at M > 6, no shock formation was observed, 
thus confirming that high Mach number shocks are sup- 
ported only when the colliding slabs have the appropriate 
temperature and density ratios. 

In conclusion, we have shown theoretically and nu- 
merically that very high Mach number shock waves are 
formed during the collision of plasma slabs. The simu- 
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lations confirmed that such shock waves, which travel 
with Mach numbers well above previous predictions 
IE S EJ E3 ' arise naturally during collision of plasma 
slabs with different electron temperatures, and driver ve- 
locity between M — 2 and M — M max . Such situa- 
tions could readily arise in astrophysics and in the in- 
teraction of high intensity lasers with plasmas ||. We 
have also shown that the maximum Mach number grows 
with the electron temperature ratio as ~ O 1 / 2 , and with 
the electron density ratio as ~ (1 + T)/T. Further 
theoretical and numerical analysis will extend the the- 
ory to extremely large Mach number shocks for which 
M > Vtcl/Vs, and include the influence of relativistic 
effects for shock velocities comparable with the speed of 
light. 
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FIG. 1: Geometry of the collisionless laminar shock wave. 
The bold line represents the electrostatic potential. The elec- 
trons from the slab L move freely, while the electrons from the 
slab R can be either free or trapped. The ions flow towards 
the shock, and are decelerated by the potential. 
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FIG. 2: Maximum Mach number of the shock as a function 
of the electron temperatures ratio O = T c r/T c l, for three 
different conditions, namely Nol = Nor (solid line), Nol = 
SNor (dash-dot line), and Nol = Nor (dashed line). The 
shaded areas represent the regions of allowed Mach number. 
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FIG. 3: Electrostatic potential ip, normalized to ksT e L/e, 
obtained from simulations (circles) and from the theoretical 
model (solid line) , as a function of the shock wave Mach num- 
ber M. 
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FIG. 4: Evolution of the ion phase space pla;l at four dif- 
ferent times steps ii = 580.8/w pe L, = 871.2/wpet, tg — 
1161.6/copeL and = 1452.0/c 1 ;p [! _t . Only the ions from the 
left slab are shown. 



